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(I) 


ABSTRACT 

The  coherent  and  incoherent  scattering  factors  and 
intensities  for  x-rays  and  electrons  have  been  calculated  for 
the  neutral  atoms  He  to  Kr  from  the  analytical  S.C.F. 
functions  of  dementi.  Results  for  x-rays  are  also  presented 
for  a  number  of  positive  and  negative  ions  of  chemical  interest. 
The  coherent  scattering  factors  can  be  used  to  assess  the  good¬ 
ness  of  these  functions  for  Evaluating  one-electron  properties. 
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I.  Introduction 


The  inherent  property  of  well-ordered  systems  like  crystals  to 
produce  diffraction  patterns  when  bombarded  with  x-ray  radiation  or 
electron  beams  has  provided  a  very  useful  method  for  the  determination 
of  chemical  structures.  In  recent  years  x-ray  analysis  has  been 
employed  with  increasing  success  in  the  structural  elucidation  of  complex 
molecules  and  macromolecules;  an  outstanding  example  of  the  latter  case 
being  the  total  structural  determination  of  vitamin  by  D.  Hodgkin 

and  co-workers  [22].  The  scattering  of  x-rays  has  also  proved  useful 
in  the  study  of  order-disorder  phenomena  and  lattice  vibrations  in 
crystals . 


Since  any  experiment  involving  the  scattering  of  electrons  or 
x-rays  necessarily  requires  the  measurement  of  intensities,  a  knowledge 
of  the  appropriate  form  factors  (German:  atomf ormafaktor )  is  necessary. 
These  factors,  more  commonly  known  as  the  scattering  factors,  can  be 
determined  if  the  electron  distribution  of  the  system  is  known. 

Tvio  methods  used  for  describing  the  electronic  distributions 
of  atoms  or  ions  are  the  Thomas -Fermi -Dirac  and  Hartree-Fock  methods. 

The  former  is  of  a  statistical  nature  and  has  been  used  to  best  advant¬ 
age  for  describing  heavy  atoms.  Corresponding  analytical  expressions 
and  values  of  the  scattering  factors  for  electrons  and  x-rays  have 
been  given  by  Bonham  and  Strand  [5], 

The  Hartree— Fock  or  self-consistent  field  method  has  also 
been  used  to  evaluate  these  quantities,  but  until  recently  only  numerical 


I 

* 

* 
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functions  of  this  type  were  available.  Calculations  of  the  scattering 
factors  using  these  functions  involved  a  tremendous  amount  of  work. 

In  i960  analytic  SCF  functions  for  helium-like  systems  were 
computed  by  Roothaan  et  al  [4l]  and  the  corresponding  coherent  form 
factors  for  x-rays  were  evaluated  by  Womack,  Silverman  and  Matsen  [50]. 

Rustgi  and  Tiwari  [43]  have  also  computed  the  coherent  x-ray  scattering 
factors  and  nuclear  magnetic  shielding  constants  for  two,  three,  and 
four -electron  systems  from  the  functions  of  Roothaan. 

Recently,  extensive  tables  of  analytic  SCF  functions  have  been 
compiled  by  Clementi  [8]  for  all  the  neutral  atoms  He  to  Kr  and  for 
a  large  number  of  positive  and  negative  ions.  These  functions  offer  an 
excellent  opportunity  to  compute  a  number  of  one -electron  properties 
like  nuclear  magnetic  shielding  constants,  quadrupole  coupling  constants# 
octopole  moments  and  so  forth.  These  quantities,  in  addition  to  others, 
have  been  calculated  for  a  large  number  of  atomic  systems  from  dementi's 
functions  [33,  34,  35] • 

It  is  the  purpose  of  this  thesis  to  give  the  necessary  formalism 
for  evaluating  the  coherent  and  incoherent  scattering  factors  and  intensities 
for  many  electron  atoms  or  ions  from  analytical  SCF  functions.  A  system- 

4 

atic  compl  lation  of  the  results  for  x-ray  and  electron  scattering  is 
presented  in  the  appendices. 

Definition  of  Terms 


The  bombardment  of  an  atom  or  an  atomic  system  by  x-ray 


. 
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radiation  will  result  in  a  scattering  of.  the  incident  wave  due  to  the 

interaction  of  the  electric  field  of  the  incident  wave  with  the  electric 

* 

field  produced  by  the  motion  of  the  atomic  electrons.  If  the  direction 

of  the  incident  wave  is  specified  by  a  unit  vector  s  and  the  direction 

o 

of  the  scattered  wave  by  a  unit  vector  s  then  a  scattering  vector  S 
is  customarily  defined  as  the  vector  difference  of  these  two? 

S  =  ?  -  s  .  (1) 

The  scattering  vector  has  a  direction  perpendicular  to  the 
plane  bisecting  the  angle  2a  between  the  incident  and  scattered  vectors. 
This  plane  is  known  as  the  reflecting  plane  since  it  reflects  the  incident 
into  the  scattered  vector,  as  can  be  readily  seen  from  the  above  definition 
or  Figure  1. 


Figure  l:  Illustration  of  the  scattering  vectors  and  scattering  plane. 


*  For  a  complete  treatment  see,  for  instance,  Kittel,  Intoduction  to 
Solid  State  Physics  [  30 3  ■> 


' 


From  the  above  figure  and  the  definition  of  the  incident  and  scattered 


vectors  as  unit  vectors  follows  the  relation 

|S  |  =  2  sin  a  (2) 

Two  scattering  sites  [e.g.  lattice  points  in  a  crystal] 
separated  by  a  distance  r  give  rise  to  a  phase  difference  [  30  ] 

5  =  (2jt/A)(r  •  S)  (3) 

between  the  corresponding  scattered  waves.  For  constructive  interference 
the  phase  difference  between  interfering  waves  must  be  an  integral 
multiple  of  2jt  . 

It  is  well-known  that  crystals  produce  diffraction  patterns 
when  bombarded  with  x  rays  of  wavelength  comparable  to  the  distances 
separating  Constituent  atoms.  Since  the  electrons  (which  are  responsible 
for  the  scattering)  are  distributed  in  a  regular  manner  about  the  atomic 
nuclei,  Debye  pointed  out  that  atoms  arranged  in  an  irregular  fashion, 
as  in  a  gas,  must  still  produce  diffraction  effects  with  x-rays.  The 
scattering  factor  is  essential  to  intensity  measurements  of  scattered 
radiation,  and  is  in  fact  a  measure  of  the  scattered  amplitude  from  a 
real  atom  relative  to  that  expected  for  a  classical  electron.  The  term 
"classical  electron"  is  taken  to  mean  a  point  charge  of  magnitude  e 
scattering  in  accordance  with  the  equation  formulated  by  J.  J„  Thomson 

[47]. 


Refer  to  James,  The  Optical  Principles  of  the  Diffraction  of  x-Rays, 
for  a  complete  discussion  [27]. 
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A  Classical  Approach  to  the  Determination  of 
Coherent  Scattering  Factors 

A  classical  method  for  determining  the  coherent  x-ray  scatter¬ 
ing  factor  will  be  considered  before  the  general  theory  is  discussed. 

An  appropriate  coordinate  system  is  arranged  as  in  Figure  2  such  that 

the  incident  and  reflected  vectors  make  equal  angles  with  the  vertical 

— > 

axis  along  which  the  vector  S  is  oriented.  This  corresponds  to  Bragg 
relection  at  an  angle  a  from  the  scattering  plane. 


Figure  2. 


The  phase  difference  between  the  radiation  scattered  by  an  element  of 


' 
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charge  located  at  (r,cp)  and  that  scattered  by  an  element  of  charge  at 

the  origin  (this  would  correspond  to  the  centre  of  the  atom)  is  given, 

«  ^  ^ 

as  before,  by  5  =  (23f/A)  (r  *  S ) .  The  ratio  of  the  amplitude  scattered 
by  a  real  one-electron  atom  to  that  scattered  by  a  classical  electron 
at  the  centre  of  an  "atom"  can  be  found  by  the  method  of  superposition,, 

If  the  charge  density  due  to  the  electron  can  be  respresented  by  a 
continuous  distribution  function  p(r),  the  scattering  factor  is  then 
given  as 

f  =  f  p(r)  e1(2*A)(?-  s)dT  . 

- >  ^ 

Since  r  describes  an  angle  cp  with  S  , 

(2jt/A)(S  •  r)  =  -^  sin  a  r  cos  cp  =  pr  cos  cp  (4) 

with  (j.  =  47t(sin  a)/A  .  If  p(r)  is  assumed  to  be  spherically  symmetric, 
then 


f  =  /  p(r )  e 


lpr  cos 


2jtr  sin  cp  drdcp 


=/• 

n 


4*r2p(r)  SiSJidr 
7 v  pr 


(5) 

(6) 


2 

Equation  [6]  is  further  simplified  by  writing  U(r)dr  =  4rtr  p(r)dr  as 
the  probability  that  the  electron  lies  in  a  spherical  shell  described 
by  radii  r  and  r  +  dr,  or 


dr 

P-r 


(7) 


■'*  ■ 
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for  the  coherent  scattering  factor.  This  factor  is  defined  as  the  ratio 
of  the  radiation  amplitude  scattered  by  the  charge  distribution  in  an 
atomic  system  to  thqt  scattered  by  a  point  electron. 

Equation  (7)  is  that  usually  found  for  the  coherent  scatter¬ 
ing  factor  in  textbooks  on  x-ray  crystallography.  Taking  the  limit  as 
(j.  approaches  zero  and  applying  L’Hospital’s  rule  to  the  integrand  of 
Equation  (7)  gives  the  familiar  result 

po° 

lim  f  =  /  lim  U(r)  5 in  dr 
“'o  u.r->0  ^'r 

~00 

=  /  U(r)dr  =  N  (8) 

J  o 

where  N  is  the  total  number  of  electrons. 

Theory  of  x-Ray  and  Electron  Scattering 

The  method  described  above  for  determining  the  scattering 
factor  is  straightforward  once  the  electron  distribution  is  known.  The 
picture  it  presents  is  misleading,  since  the  total  intensity  of  scattered 
radiation  at  moderate  scattering  angle  is  composed  of  two  parts,  coherent 
and  incoherent.  The  coherent  scattering  is  associated  with  processes 
in  which  the  scattering  system  remains  in  the  same  energy  state,  there¬ 
by  implying  that  the  scattered  and  incident  radiation  are  of  equal 
frequency.  The  incoherent  scattering  involves  transitions  from  one 
energy  state  to  another  as  a  result  of  the  perturbation,  thereby  causing 
the  scattered  radiation  to  be  of  a  longer  wavelength.  Because  of  this 
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wavelength  shift;,  the  incoherent  scattering  is  often  referred  to  as  in¬ 
elastic  or  Compton  scattering.  The  Intensities  of  both  types  of  scatter' 
ing  can  provide  useful  information  on  electron  distributions;  a  topic 
which  will  be  referred  to  the  discussion  section. 


The  total  intensity  of  x-rays  or  electrons  scattered  from  an 


atom  in  an  initial  state  described  by  ip  is  given  by 

m 


I (y)  =  V  p  Op  IP  \yj  Xip  |P *\ip  > 
v//  /  n  m1  n  1Tn  Tn  1  n  !  m 


(9) 


n 


with  the  summation  running  over  all  those  states  ip  for  which  the 


n 


excitation  energy  or  energy  difference  (E  -  E  )  is  less  than  the 

n  m 

energy  of  the  incident  beam.  The  scattering  angle  is  given  by  y  . 


For  electron  scattering  within  the  framework  of  the  first 
Born  approximation  [39] >  the  parameters  p^  are  defined  by  the  ratio 
k ^/k  ,  where  k^  and  k  are  the  amplitudes  of  the  wave  vectors  for 
the  scattered  and  incident  particles.  The  amplitude  of  the  scattered 


wave  vector  is  given  as 


k  =  [k  -  (E  -  E  )] 
n  L  v  n  nr 


1/2 


(10) 


with  the  energy  defined  per  unit  surface.  For  electron  scattering  the 
operator  is  defined  as 


p-  ['N  +  )  exp 


n 


(ii) 


k 


The  notation  of  Bonham  [4]  is  adopted.  The  operator  is  the 

complex  conjugate  of  . 


■ 
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with  the  summation  running  over  all  the  electrons  of  the  system.  The 
vector  r^  is  defined  as  the  position  vector  for  the  k'th  electron 
with  representing  the  vector  difference  between  the  incident  and 

scattered  vectors.  It  is  the  scattering  vector  similar  to  that  previously 
mentioned  for  coherent  x-ray  scattering,  but  with  ampitude 

l  l 

IS  I  =  {2k2  -  (E  -  E  )  -  2k[k2  -  (E  -  E  )f  cos  yf  (12) 

1  n 1  L  n  m  v  n  nr  ' J  v  ; 

The  classical  intensity  I  for  electrons  is  given  by  the  Rutherford 

C  Xj 

[44]  scattering  law. 

In  a  non -re lativistic  formulation,  the  parameters  p^  for 

*  .3 

x-ray  scattering  are  defined  as  vvn/v'  ’  where 

V  =  V  -  (E  -  E  )/h  (13) 

n  n  m  ' 

with  v  being  the  frequency  of  the  incident  wave. 

The  corresponding  operator  for  x-ray  scattering  is 

defined  as 


Pn  =  Yj  expfi^K^  '  ^En  "  Em^/hV^s^  '  (l4) 

j 

with  the  summation  again  extending  over  all  electrons  of  the  system. 

The  wave  number  of  an  incident  wave  (or  propagation  constant,  as  pre¬ 
ferred  by  some  authors)  with  wavelength  A  is  represented  by  k  =  2jt/A  „ 

— ) 

The  vector  S  represents  the  vector  difference  between  the  unit  vectors 

— > 

in  the  incident  and  scattered  direction.  The  vector  k  in  the  scattered 

s 

These  are  the  Breit-Dirac  recoil  factors.  For  further  information, 
refer  to  Breit  [6]. 
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direction  has  a  magnitude  which  can  be  determined  from  the  conservation 

— ) 

of  energy.  The  modulus  of  the  scattering  vector  S  is  equal  to  2  sin 
where  the  scattering  angle  7  is  again  twice  the  Bragg  angle  a  .  The 
classical  intensity  I  for  x-ray  scattering  is  given  by  the  Thomson 

C  )h 

scattering  equation  [47]. 

Since  the  coherent  scattering  only  involves  the  ground  state 
of  the  system,  Equation  (9)  reduces  to 

lc(7)  -  ^JP|V%lp*IV  -  l<VjpIV|2  h5) 

for  the  coherent  intensity.  The  operator  P  is  defined  by  Equations 
(11)  and  (14)  respectively,  with  E^  -  E^  =  0  .  By  difference,  the 
incoherent  scattered  intensity  is  given  by 

^(7)  =  1(7)  -  ic(r)  •  (16) 

The  determination  of  the  incoherent  scattered  intensity  for  electrons 
or  x-rays  usually  involves  certain  approximations  in  order  that  the 
closure  relation  may  be  used  when  summing  the  contributions  from  excited 
states.  However,  Bonham  [4]  has  developed  a  rigorous  treatment  for  using 
the  closure  relation  via  a  time -dependent  formalism  originally  suggested 
by  Lamb  [52]  and  Zemach  and  Glauber  [51,52].  The  zero-order  contribution 
given  by  Bonham  can  be  written  as 

1(7)  =  <ip  |PP*  \f  >  =  I  |P  I2  I'd,  >  (17) 

where  the  operator  P  has  the  same  definition  as  that  used  in  Equation 
(15).  Equation  (17)  depends  only  on  the  ground  state  of  the  system. 


Nc\j 
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It  represents  the  usual  Waller -Hartree  equation  for  the  total  x-ray 
scattering  or  the  Morse  expression  for  total  electron  scattering. 

From  Equations  (15,  1 6,  17)  the  incoherent  contribution  to 
the  total  intensity  is 

1.(7)  =  «p\  |p|2|^ m>  -  |<^JPI^ m>i2  •  (i8) 

1  '  m 1  1  1  1  m  1  m 1  1  m  1 

The  remaining  discussion  will  be  centered  on  Equations  (15)  and  (l8)„ 
The  notation  conventionally  employed  for  the  above  equations 


I(7)/Ici  =  f2{7)  +  S (7)  , 

(19a) 

Ic(5')/Ici  =  f2(7)  , 

(19b) 

hto/l0i  m  s(7)  ’ 

(19c) 

where  F(y)  and  S (7)  are  known  as  the  coherent  and  incoherent  scater- 
ing  factors.  They  are  defined  from  the  general  equations  above  as 

F(y)  =  <^J  ^  >  »  (20) 

j 

S (7)  =<^m||  exp{i^- r^}  |2  |#0>  -  F2(/)  .  (21) 

j 

The  theory  of  x-ray  scattering  has  been  treated  by  Waller  and 


■ 
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Hartree  [48],  Klein  [31],  and  Wentzel  [49].  A  par ticularily  useful 
equation  has  been  derived  by  Waller  and  Hartree  [48],  who  approximated 
the  total  function  ip  ,  by  a  single  Slater  determinant.  They  obtained 


— > 

expfircS 


(22) 


=  N 


X  X  (fn£& 

j  Mj 


(22a) 


with 


f 


rn}cpK(n)dT(n))5(T]y  ilk)  , 


(22b) 


where  N  denotes  the  number  of  electrons  in  the  atom  or  ion.  The 
summation  in  n  extends  over  the  electrons  while  the  summations  in  j 
and  k  run  over  the  corresponding  spin-orbitals.  The  spin  functions 
are  represented  by  q  and  • 

Equation  (22)  can  be  recast  in  the  fcsrm 


I/I  =  N  +  [  V  f..f 

'  L  4,  JJ 

j 


(23) 


since 


j 


f  f* 
jj  kk 


The  coherent  scattered  intensity  is  given  as 


,<*ss)  <  (/  P  /  [.n  Tb(rr)  q>(  *  •  th  i  ] <?xa( n j  ‘x?  \  }  .  3 


, 


-  14  - 


I  /I  =  \<lp  \  V  expfiKS* •  r*  }  hj/  >|2 
d  Q.H  1  rm 1  /  L  n .  |rm  1 

n 

(24) 

Subtracting  this  quantity  in  Equation  (23)  gives  the  Waller -Har tree 
expression 


j  j  k7j 


(25) 


for  the  incoherent  scattered  intensity.  Note  that  Equation  (25)  is 
the  expression  for  S(y)  . 

The  intensity  of  x-ray  radiation  incoherently  scattered 
increases  from  zero  with  increasing  scattering  angle,  and  approaches  the 
number  of  electrons  N  in  the  atomic  system.  On  the  other  hand,  the 
coherent  scattering  factor  (the  square  root  of  the  intensity)  decreases 
from  its  maximum  value  of  N  at  zero  scattering  angle  and  approaches 
zero  with  increasing  angle.  The  sum  of  the  two  intensities  in  an 
arbitrary  direction  is  always  constant  for  any  value  of  the  scattering 
variable,  as  shown  by  Wentzel  [491* 

The  scattering  for  electrons  is  amenable  to  a  treatment  similar 
to  that  given  for  x-rays.  In  fact,  the  electron  scattering  factors  can 
be  defined  in  terms  of  the  corresponding  factors  for  x-rays,  and  for 
this  reason  the  necessary  equations  will  be  given  in  the  next  section. 
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II.  Theoretical 


The  Formalism  for  Spherically  Symmetric 

Charge  Distributions 


The  evaluation  of  the  scattering  factors,  and  indeed  the 


* 


formalism  used  in  any  such  evaluation,  is  dependent  upon  the  model 
employed  in  describing  the  electronic  distribution  about  the  atom. 

In  the  present  treatment  recently  calculated  [8]  analytical  self- 
consistent  field  functions  were  used. 

The  scattering  vector  for  a  spherically  symmetric  atom  may 
be  taken  parallel  to  the  polar  axis  (axis  of  z-component  of  angular 
momentum)  as  the  scattering  from  such  a  system  is  independent  of  atomic 
orientation  relative  to  the  plane  of  scattering.  With  the  assumption 
of  spherical  symmetry,  the  operator  involved  in  the  calculation  of  the 
expectation  values  may  be  expanded  by  using  the  well-known  plane  wave 
expansion 


P 


for  an  electromagnetic  wave  travelling  with  its  electric  vector  parallel 
to  the  z-axis.  In  the  above  formula,  the  jp(|~ir)  are  the  spherical 
Bessel  functions,  expressible  in  terms  of  the  ordinary  Bessel  functions 
of  the  first  kind  of  order  half  an  odd  integer.  The  polar  angle  is 


designated  by  Q  while  p.  represents  the  quantity  4jt(sin  a)/A  . 


Compare  the  formalism  and  results  given  by  Bonham  and  Strand  [5]* 
who  employed  the  Thomas -Fermi-Dirac  statistical  model  in  their 
evaluation  of  the  scattering  factors. 


. 


’ 


' 

,  - 
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If  the  one  electron  state  functions  are  expanded  in  terms  of 


* 


orthonormal  basis  functions  of  appropriate  symmetry, 


<pi  - 


*tcti 


(27) 


then  a  detailed  expression  for  the  f  ^  terms  can  be  developed.  The 
following  development  is  restricted  to  the  case  when  the  basis  functions 
are  Slater-type  orbitals 


m 

xt  =  NtRt(r)Yit(e,  <P)  , 


(28a) 


with 


n  ,  r 

Rt(r)  =  r  1  e  1 


(28b) 


where  N  is  the  radial  normalization  factor,  I*  the  orbital  exponent 

m 

and  n  the  principal  quantum  number.  The  functions  Y  (0,  9)  are  the 

L  Jj 

t 

normalized  spherical  harmonics,  with  £  and  m^  being  the  quantum 
numbers  which  define  the  symmetry.  That  is,  they  are  the  azimuthal  and 
magnetic  quantum  numbers,  respectively. 

From  Equations  (26  ,  27,  28a,  28b)  it  is  readily  seen  that 

the  f.,  terms  can  be  written  in  the  form 
jk 


f.,  = 
Jk 


X  I  iP[lw(2p  +  1)]2  WtjVj 


X 


V 


For  convenience  and  simplicity  in  the  writing  of  equations,  the  symmetry 
designation  will  be  omitted. 


. 
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r°°  p2jt  o  _ 

1  Rt(r)Rv(r)j  (ur)Y  (0)Y  C(e,  <p)  X 

^  O  O  «'  O  "  p  t 


m 


Y,  (e,cp)dcp  sin  6  d6  dr  , 

JC/ 

V 


(29) 


where  the  summation  over  t  corresponds  to  cp^  while  that  over  v 
corresponds  to  cp^  •  The  above  equation  can  be  simplified  by  writing 
for  the  spherical  harmonics  the  identity 


Y>, 


1  o  X  imcp 

=  - ©(f ,  m)  e 

sl2rt 


(30) 


where  the  @($,  m)  are  the  associated  Legendre  functions  as  defined 
by  Condon  and  Shortley.  Integration  over  the  azimuthal  angle  gives 


with 


fjk  2* 


1 

X  iP(2P  +  1}  •  (2FTT]2  NtNvCtjCvk  •  <JPtv> 


p  t  V 


B(mt,  mv)  •  2jt  I  0(p,  0 )©Uv,  mt)©(iv,  n^)  sin  Q  d@  , 

o 


(31) 


r“  2 

<jptv>  =  J  r  Rt(r)Rv(r) jp(nr)dr  . 


(32) 


In  the  writing  of  Equation  (31)  the  relationship 


Y>>  -  1)2  WVT)2  e(p’  0) 


(33) 


was  employed  in  order  that  the  notation  of  the  above  mentioned  authors 
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for  integrals  over  three  associated  Legendre  functions,  viz. 

1 

d  o  ?  r ^ 

c  (i^;  i^)  =  (gp  T  i')  /  ©(p,m1-Tn2)0(i1,Tn1)0(i2,Tn2)sine  d0  (34) 

could  then  be  used  to  express  Equation  (31)  in  the  form 

fJk  =  X  X  X  iP^P  *  l)NtNvCtjCvkcP(itm’  V“)<jptv>  »  (55) 

P  t  v 

where  m  denotes  the  value  common  to  m  and  m 

t  v 

From  the  orthogonality  restriction  satisfied  by  the  eigen¬ 
functions  involving  the  azimuthal  angle  cp  follows  the  selection  rule 
that  for  a  spherically  symmetric  electron  distribution  there  is  no 
scattering  from  orbitals  with  different  z-components  of  orbital  angular 
momentum. 


Due  to  the  coupling  of  angular  momenta  the  summation  over  p 
in  Equation  (35)  Is  restricted  to  those  values  (in  steps  of  two)  for 
which  the  familiar  triangular  inequality 

|it  -  iv|  <  P  <  |it  +  iv|  (36) 

is  satisfied.  The  coefficients  of  the  radial  integrals  <j>  arising 
from  the  integration  over  6  in  Equation  (36)  are  listed  in  Table  1 
for  Slater  functions  of  s-,  p-,  d-  and  f-type  symmetry.  A  similar 
table  has  been  given  by  A.  J.  Freeman  [14],  but  it  excludes  the 
coefficients  arising  from  orbitals  of  f-type  symmetry,  which  are 
presented  here  for  completeness. 


' 
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Table 

1. 

Coefficients  of 

<jpqt> 

in 

the  matrix 

elements  f.,  for 
ik 

s,  p,  d,  and  f 

orbitals . 

Common 

Coefficient  of 

q 

t 

Coefficient 

<J0> 

<Ji 

>  <j2> 

<j^ 

<V  <j5> 

s 

s 

1 

1 

s 

P0 

is/~3 

1 

s 

do 

/5 

-1 

s 

£o 

isF] 

-1 

po 

po 

1 

1 

-2 

po 

do 

3i/l5 

2 

-3 

po 

fo 

^3/7 

-3 

4 

do 

do 

1 

1 

-10/7 

18/7 

do 

£o 

i/3^35 

27 

-28 

50 

fo 

fo 

1/33 

33 

-44 

54  . 

-100 

p-.i 

P+1 

1 

1 

1 

p+i 

dn 

31//5 

1 

1 

p+i 

f+l 

5/177 

-1 

-1 

d+i 

d+l 

1 

1 

-5/7 

-12/7 

d+i 

f+l 

18 

-7 

-25 

f+l 

f+l 

1 

1 

-1 

3/11 

25/U 

d+2 

d+2 

1 

1 

10/7 

3/7 

d+2 

f+2 

i/5/7 

9 

14 

5 

f+2 

£-* 

1 

1 

0 

-21/11 

-10/11 

f+3 

f+3 

1 

1 

5/3 

9/n 

5/33 

,  h  ,  q 
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Methods  for  Treating  Noa-Spherical  Charge  Distributions 

If  the  electron  distribution  is  no  longer  spherically  symmetric 

(e.g.  open  shells  of  electrons  of  p-  or  d-type  symmetry),  the 

scattering  will  be  dependent  upon  the  orientation  of  the  scattering 

vector.  Equation  (35)  is  no  longer  strictly  valid,  as  the  f  terms 

— > 

should  be  evaluated  for  an  arbitrary  orientation  of  S  . 


Methods  for  handling  this  problem  are  discussed  in  the 
literature.  McWeeny  [37 3  has  demonstrated  that  the  scattering  from  an 
orbital  of  p-type  symmetry  oriented  in  an  arbitrary  direction  can  be 
completely  described  by  two  factors  f||  and  f_j_  ,  and  an  angle  (3 
between  the  orbital  axis  and  the  vector  S  .  In  this  notation  fj_ 
represents  the  scattering  factor  for  an  orbital  of  p-type  symmetry  with 


axis  oriented  perpendicular  to  S  ,  while  f||  refers  to  the  scatter- 

— > 

ing  factor  when  S  is  parallel  to  the  orbital  axis.  Due  to  geometric 

considerations,  this  method  is  not  easily  generalized  to  orbitals  of 

higher  angular  momentum.  Furthermore,  if  an  average  is  carried  out  over 

— > 

all  orientations  of  S  ,  e.g.,  an  integration  over  the  angle  (3 
weighted  by  an  element  of  solid  angle,  2jt  sin  (3  d(3  ,  the  result  is 
equivalent  to  that  expected  on  the  assumption  of  spherical  symmetry. 

In  other  words,  the  scattering  from  an  open  shell  of  p-electrons  is 
proportional  to  that  expected  from  an  equivalent  closed  shell  of  p- 
electrons,  the  constant  of  proportionality  being  the  ratio  of  the 
occupation  numbers. 


More  recently,  Freeman  [14]  has  considered  the  evaluation  of 
the  scattering  factors  for  open  shells  of  p-  and  d-type  symmetry  by 


1 
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group  theoretical  methods.  An  orbital  with  its  axis  oriented  in  an 
arbitrary  direction  can  be  brought  parallel  to  the  z-axis  by  a  coordinate 
transformation  of  the  basis  functions.  The  transformation  involved  is 


a  rotation, 


through  the  three  Euler  angles.  As  the  radial 


parts  of  the  basis  functions  are  invariant  under  the  group  of  three- 
dimensional  rotations,  only  the  rotation  of  the  spherical  harmonics 
themselves  must  be  considered,  a  task  complicated  by  the  fact  that  the 
addition  theorem  for  spherical  harmonics  is  no  longer  satisfied.  The 
results  derived  by  Freeman  for  p-type  orbitals  are  analagous  to  those 
given  by  McWeeny;  in  addition,  Freeman  lists  formulae  for  matrix  elements 
involving  orbitals  of  d-type  symmetry. 

Freeman's  method  clearly  exhibits  the  angular  dependence  of 
the  matrix  elements  f  ^  on  the  orientation  of  the  scattering  vector, 
but  it  has  the  disadvantage  in  that  the  orbitals  in  the  rotated  system 
are  no  longer  eigenfunctions  of  the  orbital  angular  momentum  operators. 
This  follows  from  the  fact  that  the  unitary  rotation  operator  R  does 
not  commute  with  these  operators.  Hence  they  cannot  share  a  set  of 
eigenfunctions  with  R  . 

The  consequence  of  this  non-commutation  can  be  demonstrated 
in  the  following  manner  .  If  G  is  a  unitary  operator  corresponding 


to  the  transformation 


and  |>  is  an  arbitrary  ket  from  the  domain 


of  G  ,  then 


(37) 


* 


See,  for  instance,  Messiah,  Quantum  Mechanics,  Vol.  2,  page  1089° 


' 
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and  if  0  is  an  observable  of  a  quantum  system  then  in  the  transformed 
system 


Jrm  - 


GOG 


-1 


(38) 


Consequently,  if  L  =  (Lx,  L^,  L^)  is  a  vector  operator  attached  to 


f  <  7 

the  system,  L.  =  L  *e.  ,  with  e.  a  basis  vector 


then 


A 


1  u 

[L.  ]  =  GB.G_i  =  (L  •  A.  )  =  G,  .B  . 
i  1  1  ij  J 


(39) 


where  G^  represents  the  matrix  element.  By  designating  L  as  the 
momentum  operator  and  as  the  rotation  in  the  above  discussion,  it 

follows  that  the  projections  of  the  orbital  angular  momentum  in  the 
rotated  system  are  modified  by  the  matrix  elements  of  the  rotation 
matrix. 


Averaging  Procedure  for  Aspherical  Charge  Distributions 


The  method  to  be  described  here  for  the  evaluation  of  the 
scattering  factors  for  non-spherical  charge  distributions  does  not 
exhibit  any  angular  dependence  since  an  averaging  procedure  originally 
outlined  by  Pohler  and  Hanson  [40]  is  used.  Instead  of  employing  the 
plane  wave  expansion  for  a  wave  travelling  along  the  positive  z-axis, 
consider  the  expansion 


00  +p 

exp{i(j.r  cos  0}  =  4jt 

p=0  rnt-p 


ipjp(nr)Y"(e,(p)Y^n(e',  r) 


(40) 


* 


A.  R.  Edmonds  Angular  Momentum  in  Quantum  Mechanics  [12]. 


■ 
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for  a  wave  travelling  in  an  arbitrary  direction  specified  by  the  coordinates 

f  *  ml  * 

6,  cp  .  The  average  incoherent  scattering  factor  for  an  aspherical 
electronic  charge  distribution  is  determined  by  integrating  (averaging) 
the  expression  for  the  f  terms  over  the  surface  of  a  sphere,  viz. 


4tc 


f  1^  sin  6  * 
Jk  1 


deT  dcp T  =  < 


(41) 


For  orbitals  expanded  in  terms  of  Slater-type  functions,  the 
maxtrix  elements  are  now  defined  by 


(42) 


Introduction  of  the  identities  for  the  spherical  harmonics 
into  Equation  (42)  and  integration  over  the  azimuthal  angle  gives 


1 


p  t  v 


* 

m  -  m 

x  <jptv>  Yp  V(e,,cp’)  (43) 


for  the  individual  matrix  elements.  The  notation  for  the  radial  and 


-x- 

In  what  follows,  the  term  aspherical  is  taken  to  mean  non-spherical « 


, 
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angular  integrals  is  the  same  as  that  employed  before.  Multiplication 
of  Equation  (43)  by  its  complex  conjugate  and  integration  over  all 
angles  Qx ,  Cp’  f  making  use  of  the  orthonormality  conditions  satisifed 
by  the  spherical  harmonics,  leads  to  the  expression 


<  f 


jk 


£  ,2  . 
f  sin  6 

jk  1 


de'  dcp 1 


I 

p  t  v 

for  the  average  of  the  square  of  the  individual  matrix  elements.  The 
summation  over  the  index  p  is  again  restricted  to  those  values  which 
satisfy  the  triangular  relationship.  Inspection  of  Equations  (44)  and 
(35)  indicates  an  important  difference  concerning  the  parameters 

p 

C  (ijm^.;i^m^)  .  In  Equation  (35)  only  those  terms  remain  for  which 
mt  =  m^  ;  a  selection  rule  which  no  longer  holds  in  Equation  (44). 

The  interpretation  can  be  given  in  the  following  manner.  By 

using  the  plane  wave  expansion  and  integrating  over  all  orientations  of 
— ) 

S  ,  the  region  of  space  over  which  the  orbital  was  originally  defined 
has  been  generalized  to  a  sphere.  Since  this  has  the  physical  effect 
of  "smearing"  the  corresponding  electron  distribution  throughout  a 
spherical  volume,  the  averaging  procedure  is  equivalent  to  an 
assumption  of  spherical  symmetry. 

The  averaging  method  should  give  the  same  results  as  the  one 


,2„2„2  „2 


(2p+l)NtNvCt.CvkC 


2  2 
vmv)  <jptv> 


(44) 


derived  on  the  assumption  of  spherical  symmetry.  The  obvious  way  to 


.  - 


. 
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check  the  equivalence  of  the  two  methods  would  be  to  determine  the 
coef ficients  of  the  radial  integrals  for  an  arbitrary  f  ^  term. 

This  is  readily  accomplished  by  using  the  selection  rule  and  the  square 
of  the  absolute  value  of  the  appropriate  coefficients  listed  in  Table  1 
and  comparing  the  results  obtained  to  those  found  by  using  the  numerical 

p 

values  for  the  C  (i^ m ^ °f  Equation  (kk)  weighted  by  the  values 
of  (2p  +1)  .  Naturally,  the  values  found  by  both  methods  is  dependent 
upon  the  occupation  number  of  the  shell. 


In  connection  with  the  scattering  from  open  shells,  only 
those  terms  which  contribute  to  the  incoherent  scattering,  i.e»,  the 
f  ^  >  have  been  considered.  If  an  average  over  aspherical  charge  dis¬ 
tributions  or  spherical  symmetry  is  assumed,  the  terms  contributing  to 
the  coherent  scattering,  the  f  ,  will  contain  radial  integrals 
involving  only  the  zero-order  Bessel  function  regardless  of  the  symmetry 
or  occupation,  number  of  the  shell.  The  coefficients  of  higher  order 
terms  always  cancel,  as  can  be  easily  verified  from  the  coefficients 
listed  in  Table  1.  To  illustrate  this  point,  consider  the  contribution 
to  the  coherent  scattering  from  a  d-shell  containing  n  electrons, 
n  <  10  . 


Since  an  average  over  an  aspherical  electron  distribution  or 
a  spherical  distribution  has  been  assumed,  the  scattering  from  this  shell 
will  be  directly  proportional  to  that  expected  from  a  closed  d-shell. 


=  ^{2fo(d)  +  4f1(d)  +  Uf2(d)} 

=  15  {2(i-<j0>  -  Y  '<J2>  +  T  +  1*(1'<Jo> 


E»«  aria  gel  >oo3*  ^Jlbaoa  a  1  alriT 


■ 
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+  4(i-<jo>  +  Y  -<j2>  +  2  <ju>)) 

=  T5  (I0<j0»  (45) 

=  n<jQ>  ,  n  <  10  . 

It  is  this  absence  of  higher  order  and  cross  terms  (terms  involving 
orbitals  of  different  symmetry  or  principal  quantum  number)  which 
facilitates  a  relatively  simple  evaluation  of  the  coherent  factors  and 
intensities . 


Evaluation  of  the  Radial  Integrals 

The  evaluation  of  the  radial  matrix  elements  will  now  be 
considered.  The  n-th  order  spherical  Bessel  function  may  be  reduced  to 
a  polynomial  in  r,  involving  sin^s  and  cosines.  By  a  modification  of 
a  recursion  formula  given  by  Hildebrand  [21], 

ji>r>  *  Jn  -  A0  '  jn  -  2(l"r)  (U6a) 

and  the  two  basic  relations 

j0(n*)  =  (sin  |ar)/|ir  (46b) 

j_1(|ir)  =  (cos  p,r )/p.r  (46c) 

the  integral  <jpqt>  can  be  written  as  a  combination  of  integrals  of 


the  general  form 


’ 
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I  =  /  R  (r)R  (r){Sin  ^r}  (— )r2dr 

Jo  tK  '•cos  |arJ 


r°°  "t  +  "v  -(^t  +  ^V)r  rsin  ^r)  (— )dr 
-lx  e  '•cos  )j,rJ  'jar7 

n 


These  integrals,  written  in  general  as 


*■*>  -  i  /: 


n  -er  , 

r  e  sm  jj.rdr  , 


_00 

o 


n  -er  , 

r  e  cos  [irdr  , 


can  be  readily  evaluated  by  recognizing  that 


P( 

S>r)  =  (-l)n  / 

0  o 


oo  n 

O  r  -er  •  i, 

Le  sm  |-irjdr 


be 


n 


=  (-D 


n  3 


n  o,oo 


■n  n 

be  J  o 


-er  , 

e  sin  |ar  dr 


n 


N  n  2  2 

e  +  |j. 


] 


and 


°n  +  Ar)  =  nfcr) 


Repeated  differentiation  of  Equation  (49)  and  collection  of  the 


* 


yields  the  general  form  , 


See  also  the  work  of  Rustgi  and  Tiwari  [43]  • 


(47) 


(48a) 

(48b) 


(49) 


(50) 

terms 


■ 


with 


C  n  +  1 
2  ’ 


n*  =  v 


n  i 

2  +  1  ’ 


for  n  odd 


for  n  even 


Carrying  out  the  operation  indicated  by  Equation  (50)  then 


gives 


n1 


Cn+l(^r) 


y  (-i)p+i  (2P1)'*2(p'i)  en'2(p"i} 

* }  P=i 


(c2+|l‘) 


2  2 

{A|_i  +  Be  J 


with 


A  =  2(p  -  n)  -  3  , 

B  =  2p  -  1  , 

nr  taking  the  values  previously  defined.  The  formal  interchange  of 
differentiation  and  integration  in  Equation  (49)  is  justified,  since  the 
integrand  is  continuous  and  defined  everywhere  on  the  interval  (0,°°)  „ 


Sample  Formulae  For  Krypton 


As  an  example  of  the  application  of  the  formulae  developed  thus 
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* 


far,  consider  the  evaluation  of  the  coherent  and  incoherent  scattering 
factors  for  Krypton.  The  detailed  expression  for  the  summation  terms 
in  Equation  (25)  is  given  by 


ore2  c2  e2  e2  1  ar2  e2  e2  1 

2  [  f  +  f  +  f  +  f.  ]  +  6[f0  +  f^  +  f.  ]  + 

Is  2s  3s  4s  2p  3P  4p 

1ne2  ,  r-2  -2  J2  2  .  ^2 

3d  +  Is, 2s  +  Is,  3s  4  Is,  4s  +  2s,  3s  +  2s, 4s 


'd  1  Tore2  e2  e2  e2  e2 

'f  f7  \  J  4*  12  If-  0  -f-f-  f  -  1  -{*  ^ 

3s, 4s  ls,2p  ls,3p  ls,4p  2s, 2p  2s, 


3p 


2  2  2  2  2  2 

2s, 4p  +  3s, 2p  '  3s, 3p  +  3s, 4p  +  4s, 2p  *r  4s, 3p 

e2  1  onl-e2  e2  e2  e2  1 

f,  |  J  +  201  f.  H*  f  ~  +  f-.  J  4"  f|  J  J  4" 

4s,4p  Is, 3d  2s, 3d  3s, 3^  4s, 3d 

2  2  2  2  2 

12  [ff  0  +  fZ  _  +  ff  ,  +  ff  ,  .  +  ff  f  r\  \  I  + 

2p,2p  3P,3P  4p,4p  2p(0)3p  2p(0)4p 


(53) 


2  2  2  2 
f3p(o)4p^  *  2^^f2p(2)3p  ^  f2p(2)4p  4  f3p(2)4p^  + 

24[f2p(l)3d  +  f3p(l)3d  +  f4P(l)3d]  +  36[f2p(3)3d  + 

2  2  2  ,  100  2  180  2 
3p(3)3d  +  3p(3)3d  +  4P(3)3dJ  +  7  3d(2)3d  +  7  3d(4)3d 


Reference  must  be  made  to  the  calculations  of  James  and  Brindley  [28] 
by  an  interpolation  procedure.  Calculations  have  also  been  carried 
out  by  Berghius  et  al  [3],  Keating  and  Vineyard  [29],  Hurst,  Miller 
and  Matsen  [24],  Freeman  [15,  16],  Silverman,  Platas  and  Matsen  [45], 
Hurst  [23],  Womack,  Silverman  and  Matsen  [50],  Bonham  and  Strand  [5], 
Rustgi  and  Tiwari  [43],  Cornille  and  Roux  [11],  and  Polher  and  Hansen 
[40],  among  others. 


The  writing  of  this  and  other  expressions  is  simplified  by  the 
grouping  of  similar  terms  with  the  common  factor  explicitly  written. 
The  notation  for  the  scattering  factors  is  defined  below.  Subindices 
denote  orbitals. 
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The  number  of  terms  contributing  to  the  incoherent  scattering  factor  for 
a  system  with  s-,  p-,  and  d-type  electrons  is 

\  n(n  +  1)  np(np  +  l)  +  npnd  +  2nd  (54) 

where  n  is  the  total  number  of  occupied  orbitals  and  n  and  n, 

p  d 

are  the  number  of  occupied  orbitals  of  p-  and  d-type  symmetry, 
respectively.  The  coherent  scattering  factor  is  given  by  Equation  (24), 
as 

F  -  £[£ls  +  £2s  +  f?s  +  f4s]  +  6t£2p  +  f3p  +  V  +  l0£Jd  (55) 

Notation  for  the  Scattering  Factors 

The  notation  employed  in  Equations  (53)  and  (55)  is  that 
originally  used  by  Keating  and  Vineyard  [23],  extended  later  by  Freeman 
[14],  with 

£  -  y 

ns 

t  v 

f  ,  =  V  V  N«-N  c  C  /  R  (r)R'(r)j  (pr)r2dr  , 

ns,n's  /  tv  tns  vo.' s  J  tv  '  vv  /JCr^  '  ’ 

t  v 

f  =  V  V  N  N  C  C  /  R  (r )R  (r ) jA(|ir )r2dr  , 

np  Z_i  Z_i  v  tnp  vnp  J  t  v  0^  '  ’ 

t  v 

f  =  )  )  N  N  C  C  /  R  (r )R  (r) j_(^r)r2dr  , 

np,np  tv  tnp  vnp  J  t  v  J2vr/ 

t  v 


N  N  C  C 
t  v  tns  vns 


J "  Rt(^)Rv(r)j0(M-r)r2dr  , 
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'np^jn1  p 


)  )  N  N  C  C  /  R  (r)R'(r)j  (|j,r)r2dr  ,  m  =  0,2 

/  /  tv  tap  vn'p  J  t'  '  v  Jnrr  *  ’ 


t  v 


ns,n'p 


V  V  NNC  C  ,  r  Rfr^'O^^r^dr 

tv  tns  vn'p  J  tv  '  v  '  ’ 


t  v 


'np(m) 


,  ,  =  y  y  N  N  C  C  ,  ,  /  R  (r )R ' (r ) j  (^ir)r2dr  ,  m  =  1,3 
n’d  Z_i  Z_i  t  v  tnp  vn'd  J  t  v  Jm^  y 


t  v 


nd 


X  1  NtNvCtndCvnd  /  Mr)Rv(r)j>r)r2dr 


t  V 


‘nd  (m)n'  d 


X  X  Ct»dCvn'd  f  Rt(r)^(r)jm(Mr)r2dr  ,  m  =  0,2,4 


t  V 


(56) 


The  summations  in  the  above  expressions  extend  over  the  corresponding 
basis  functions.  The  order  of  the  spherical  Bessel  functions  defines 
the  individual  factors  in  an  unambiguous  fashion. 


A  Sample  Calculation 

In  order  to  illustrate  the  evaluation  of  one  of  the  above  terms, 

consider 

£3d(4)5d  =  X  X  VvCt3dCv3d  /  Rt(r)Rv(r)V^r)r2dr 

t  V 

Y,  NtNvCt3dCv3d  *  f105  /  rne"er(sin  ^r)/([ar)5dr 
t  v 


I  ■  bn* 
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-45  J'  rne  £r(sin  jj.r  )/(p.r  )^dr  +  J' rn  e  er(sin  (j.r)/(ar  dr 

-  105  J'  rne  €r(cOs  )ar )/(iar )^dr  +  10  J' rn  e  £r(cos  p.r  )/(p.r  )2dr  ] 

-  3840  £  £  Wt3dCv3df^/(^  +  42)6)  • 

t  v 

For  this  case,  n  =  n  +  n  =6  for  all  t,  v  since  the  3-d  orbital 
for  Krypton  was  expanded  in  terms  of  Slater  orbitals  of  this  symmetry 
and  quantum  number.  The  sum  of  the  orbital  exponents,  £  +  £v  ,  has 
been  designated  by  the  variable  e 

Analytical  expressions  such  as  the  one  given  above  were 
incorporated  in  a  general  program  and  evaluated  by  computer.  Values  are 
presented  in  the  appendix  for  (He  to  Kr)  for  the  incoherent  atomic 

scattering  factors  S(|j.),  the  coherent  atomic  scattering  factor  F((jJ, 

2  2 

the  coherent  scattered  intensity  F  (p.)  and  the  quantity  [N  -  F (|_l )  ] 

for  a  wide  range  of  values  of  the  scattering  variable,  |_i  =  4jt(sin  <x)/A  , 

in  1  „  Besides  the  individual  factors  f .  .  which  necessarily  had 

1  j 

+ 

to  be  evaluated,  values  of  the  sums  of  the  squares  of  the  diagonal 
and  off-diagonal  elements  and  the  ratio  of  the  incoherent  to  coherent 
x-ray  intensity  have  been  evaluated  for  all  the  neutral  atoms  He  to 
Kr  and  for  a  number  of  positive  and  negative  ions  of  chemical  interest. 


(57) 


* 


Unpublished  results. 
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Scattering  Factors  For  Electrons 


The  electron  scattering  has  been  treated  in  a  manner  similar 


to  that  described  for  x-rays.  The  operator  yielding  the  expectation 

•X--X- 

values  was  given  as 


Q  =  (-N 


exp  { i  kS 


r  }) 
nJ 


(58) 


with  the  classical  intensity  defined  in  a  manner  different  from  that 
for  the  x-ray  scattering.  Making  use  of  the  notation  introduced  above 
for  the  coherent  and  incoherent  x-ray  scattering  factors,  the  total 
scattered  intensity  for  electrons  can  be  written  as 


I/I  =  N2  -  2NF  +  F2  +  S 
/  eg, 

=  (N  -  F)2  +  S 


(59) 


2 

with  the  coherent  intensity  defined  as  I  /I  =  (N  -  F )  .  The  coherent 

c  C  g 

scattering  factor,  customarily  designated  as  f^(|j.)  >  is  then  the  square 
root  of  this  quantity  except  for  the  case  when  the  scattering  variable 
p,  assumes  the  value  zero.  By  a  limiting  procedure,  the  expression 
known  as  Ibers  [25]  formula 

f/o)  =  i  <r2>  (60) 


can  be  obtained.  This  equation  has  been  rederived  by  Strand  and  Bonham 


For  details,  reference  can  be  made  to  the  work  of  Bartel  and  Gavin,  Jr, 
[1,2],  and  to  that  of  Iijima,  Bonham  and  Ando  [26], 

Equation  [11]  of  the  theory  section. 


1 
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[46].  Table  1  of  Appendix  I  lists  the  values  of  f^^(0)  for  a 
series  of  positive  ions,  neutral  atoms,  and  negative  ions.  These 

fg  (o)  values  were  determined  from  Equation  (59)  using  the  total 

2 

expectation  value  of  <r  >  determined  from  the  analy  -fcical  Hartree- 
Fock  functions  of  Clementi  [8]  and  Malli  [36].  Values  for  the  neutral 
atoms  only  have  been  given  by  Bonham  and  Strand  [5]  and  Strand  and 
Bonham  [46] . 
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III*  Discussion 

The  averaging  procedure  described  above  is  valid  for  the 
scattering  from  a  system  of  non-interacting  gaseous  atoms  whose  space 
axes  are  randomly  oriented  with  respect  to  the  scattering  plane »  The 

assumption  of  spherical  symmetry  does  not  conflict  with  the  picture 

■jf 

presented  by  authors  of  reference  works  on  x-ray  diffraction^  wherein 
the  x-ray  is  described  as  interacting  with  the  time-average  electron 
cloud  of  the  scattering  system  and  setting  it  into  sympathetic  motion,, 

In  the  absence  of  any  forces  tending  to  remove  the  degeneracy  (e„g.; 
Zeeman  or  Stark  effects),  the  time-average  distribution  of  a  degenerate 
open  shell  will  be  spherically  symmetric. 

The  subject  of  electron  and  x-ray  scattering  is  directly 
related  to  the  general  problem  of  scattering  by  a  central  potential. 
Equations  for  electron  and  x-ray  scattering  are  an  outgrowth  of  this 
more  general  case,  and  as  such,  are  subject  to  the  same  approximations 
involved  in  the  general  problem.  For  instance,  the  summation  in 
Equation  (9)  for  the  electron  scattering  involved  only  those  excited 
states  for  which  the  energy  of  the  incident  beam  was  smaller  than  the 
energy  differences,  En  -  Em  ,  between  the  ground  and  excited  states. 
This  is  nothing  more  than  a  manifestation  of  the  first  Born  approximation, 
but  it  insures  that  the  incident  electron  will  not  become  bound.  The 
problem  of  x-ray  scattering  involves  the  calculation  of  the  amplitude 
of  scattered  radiation  when  an  electromagnetic  wave  interacts  with  a 

* 


See,  for  example,  James,  The  Optical  Principles  of  the  Diffraction 
of  X-rays  [27], 
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system  whose  unperturbed  wave  functions  are  known.  To  deal  satisfactorily 

with  this  problem,  the  total  radiation  field  and  scattering  system 

should  be  considered  as  forming  a  single  quantized  system,  as  proposed 
* 

by  Dirac.  Considerations  such  as  these  are  beyond  the  scope  of  this 
work.  Rather,  some  interesting  applications  of  x-ray  scattering  data 
will  be  examined. 

Determination  of  Distribution  Functions  From  Scattering  Data 

Conventional  x-ray  crystallography,  based  on  measurements  of 
coherently  scattered  radiation,  is  a  method  used  for  determining  chemical 
structures.  However,  the  measurement  of  coherent  scattered  intensities 
can  be  applied  to  another  problem,  namely  the  determination  of  radial 
charge  density  distributions.  The  coherent  scattering  factor  is  given 
as  the  expectation  value  of  a  one-electron  operator  and  can  be  used  to 
describe  the  electron  density  about  nuclei,  a  one-electron  property. 

This  is  readily  seen  from  Equation  (15)  which  reduces  to  the  familiar 
expression 

Ic(n)/Ici!  -  f2(h)  ,  (61) 

with 


for  spherically  symmetric  atoms  or  an  average  over  random  orientations 
of  aspherical  atoms,  as  previously  shown.  Again,  D(r)  is  the  radial 
distribution,  4jtr  p(r),  of  planetary  electrons,  with  p(r)  the  mean 

* 


Reference  can  be  made  to  Fermi  [13l» 


■ 
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electron  density  function  with  respect  to  the  radial  distance  r  from 
the  nucleus.  In  Equation  (62)  the  function  D(r)  can  be  regarded  as 
as  sum,  ^  D^(r)  ,  of  individual  one-electron  distributions. 

T 

Equation  (62)  very  clearly  indicates  that  the  electron-nuclear 
radial  distribution  D(r)  is  related  to  the  coherent  scattering  factor 
(as  determined  from  experimental  intensity  measurements  of  coherently 
scattered  radiation)  by  a  Fourier  sine  integral.  Taking  the  indicated 
transform  gives 

p  r00 

D(r)/r  =  -  /  |i  F(|_i)  sin  jar  d|i  (63) 

J  o 

with  | a  =  4it(sin  a)/A  .  The  differentiation  between  the  elastic  and 
total  intensity  is  rarely  carried  out  experimentally,  but  it  can  be  done 
in  principle.  Compton  [9]  has  accomplished  this  differentiation  for  the 
larger  scattering  angles. 

■£ 

Information  on  the  elctron-electron  radial  distribution,  a 

two-electron  property,  can  be  determined  from  the  total  scattered  intensity. 

That  the  operator  giving  the  total  scattered  intensity  depends  on  the 

coordinates  of  two  electrons  is  recognized  by  writing  the  Waller  Hartree 
** 

equation 


— 7 

exp{ijrcs 


?  }f\f  > 

nJ  1  1  m 


(64) 


Refer  to  Bartell  and  Gavin  [1,2]  for  further  details. 
Equation  (22)  of  the  theory  section. 


■ 


. 
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for  the  total  intensity  in  the  form 


--  <U  Y,  Y  [expfiK®*‘  ^HexpfiK?. ^k)]*l^m> 


<^1  y  exp(iKs>-(?j  -\)}i^m> 


(65) 


<U  Y,  X  ' 


with  t .,  s  r .  -  r  .  Rather  than  taking  the  3N  spatial  coordinates 
J  k 

for  the  N-electron  wave  function  ip  from  the  components  of  the  vector 

m 

— >  — >  — >  — >  — > 

set  [r  ,  .  .  .  ,r  .  ,r  ,r  ,  ...,r  J,  the  components  of  the  set 
I  J  k  £  h 

— >  — >  — >  ' — > 

[r  ,  . . .  ,r  . , t  ,r  ,  . . .  ,r  ]  can  be  employed.  The  wavefunction  so 
1  J  J  k  &  n 

constructed  then  implies  the  electron-electron  distribution  function, 

P  (t..  )  for  the  i'th  and  k’th  electrons.  It  can  be  defined  by 
n\  j  j 


-  f  KK  dT 


(66) 


where  the  integration  is  to  be  carried  out  over  all  coordinates  except 


those  of  t 


jk  ’ 


Physically,  Pn(tjk)dTjk  (dTJk  =  t jkdt jfc  sir,  0jk  d6jk  d9jk) 

can  be  interpreted  as  a  probability.  It  describes  the  probability  that 

the  distance  separating  electrons  j  and  k  lies  between  t  and 

t  ^  +  <^tjk  t*ie  P°lar  anc*  azimuthal  angles  defining  the  orienta- 

*— > 

tion  of  the  vector  t lie  between  Q and.  6.,  +  d 6.,  and  between 

JK  jK  JK  JK 


q» 


■ 

' 
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Cf>jk  and  +  dtp^  ’  respectively.  The  angles  are  reckoned  from 
any  conveniently  chosen  reference  frame.  These  arguments  imply  that 
the  total  scattered  intensity  can  be  expressed  in  the  form 

Vhz  =  X  X  /  e*P{1“-t>jk3dTjk  ’  (67> 


which  simplifies  to 


V1 


c  Z 


f  pjk(t)  sln  dt 


(67a) 


when  the  atom  is  spherically  symmetric  or  when  an  average  over  random 

orientations  of  s  is  considered.  In  Equation  (67a),  P.,(t)  is  the 

Jk 

radial  distribution  function  of  electron  j  with  respect  to  electron 

k  .  The  mean  density,  p.,  (t)  of  one  electron  with  respect  to  another 

Jk 

is  defined  by  the  delation 


Iwr  Pjk(t)  =  Pjk(t) 


(68) 


and  a  total  electron-electron  distribution  function  by  the  relation 


P 


V(t) 


(69) 


where  5(t)  is  the  usual  Dirac  delta  function  and  N  is  the  number  of 
electrons.  The  expression  for  the  total  scattered  intensity  is  now 


■ 
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given  as 


VIci  -  t  f  P(t)  sin  dt  +  N  }  .  (70) 

Taking  the  Fourier  transform  of  Equation  (70)  then  gives 

pO° 

p(t)/t  =  f  /  VXC  »  -  sin  tit  dp.  (71) 

4  o 

for  the  total  electron-electron  distribution  function.  Lack  of  experi¬ 
mental  intensity  measurements  to  (sin  <x)/A  =  00  f  which  are  required 
in  both  Equations  (63)  and  (71),  can  be  treated  by  a  procedure  suggested 
by  Kauptman  and  Karle  [20] . 

Equation  (71)  suggests  a  method  for  determining  a  two- 
electron  property  experimentally.  The  function  P(t)  is  difficult  to 
calculate  theoretically  and  it  seems  reasonable  that  an  experimental 
measure  of  the  total  intensity  would  afford  a  more  accurate  deter¬ 
mination  of  P(t)  than  would  a  calculation  based  on  one-electron 
wavefunct ions . 

Determination  of  Electronic  Energies  from  Scattering  Data 

Total  electronic  energies  can  be  found  from  the  measurement 
of  scattered  x-ray  intensities,  but  in  the  light  of  better  alternative 
experimental  methods  the  results  may  not  be  competitive.  The  method 
is  described  in  the  following  manner.  The  non-relativistic  Hamiltonian 
operator  for  an  atom  or  ion  is 


- 


Sg-g.i 
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H  =  A  [t  +  Vj  +  A  A  V 

j  j  <  k 


(72) 


where  TA  and  U  represent  the  kinetic  energy  and  electron-nuclear 
potential  operators  with  U  ^  the  electron-electron  repulsion  operator 
Thus  the  average  potential  energies  are  simply 


*00 


U  .  =  -e 
PJ 


[D  (r)/r]dr 


(73) 


=  -e 


J '  [D(r) j/r  ]dr 


and 


j  MJ 


U  =  e 
jk 


r  1 1 .  V** 


j  k/J 


(7*0 


=  e 


r°° 

/  [p(t)/t]dt 
o 


with  the  total  energy  given  by  the  virial  theorem, 


E  =  £  O  U, 


2  L  L  p1  L 

P  j  W 


U.,  ) 
jkJ 


(75) 


d  p°0  pOO 

=  \  {  [p(t)/t]dt  -  /  [D(r)/r]dr) 

^  o  o 


The  total  energy  calculated  in  this  fashion  accounts  for  the  correlation 
energy,  which  is  difficult  to  determine  theoretically. 


; . 


Validity  of  Results 


The  accuracy  of  the  expectation  values  calculated  will  now  be 
considered.  In  order  not  to  become  involved  in  a  tabulation  of  values, 
either  theoretical  or  experimental,  the  discussion  of  the  present  results 
will  mainly  depend  on  a  theoretical  rather  than  a  comparison  approach. 

As  pointed  out  by  Goodisman  and  Klemperer  [18],  an  approximate 
wavefunction  0  used  in'  describing  a  system  always  gives  the  energy  to 
a  higher  degree  of  accuracy  than  other  properties,  since  the  true  wave- 
function  Y  is  an  eigenstate  of  the  energy.  The  true  state  function 
¥  can  be  written  in  terms  of  the  approximate  function  $  and  a 
correction  function  x.  chosen  orthogonal  on  the  space  of  $  , 

<X*|$>  =  0  ,  by  setting 

2 

*  =  (1  -  +  <=x  •  (76) 

The  parameter  e  is  a  measure  of  the  deviation  of  the  approximate 
solution  from  the  exact  solution,  and  always  takes  the  values  0  <  €  <  1  . 
It  rapidly  approaches  zero  for  good  approximate  functions  0  . 

Any  one-electron  property  calculated  by  an  approximate  function 

might  then  be  expected  to  contain  a  relative  error  of  magnitude  e  , 

except  for  the  Hartree-Fock  functions  where  the  error  is  proportional 
2 

to  e  .  This  is  a  consequence  of  Brillouin’s  [7]  theorem  which  states 
that  the  first-order  correction  to  the  single  determinant  Hartree-Fock 
function  describing  a  closed  shell  system  is  orthogonal  to  the  Hartree- 
Fock  function  in  the  coordinates  of  two  electrons.  This  means  that 


' 


-!3  , 
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any  correction  to  the  expectation  value  of  a  one-electron  operator  Q 

computed  from  a  Hartree-Fock  function  $  will  only  contain  terms  in 

2 

£  and  higher.  The  correction  to  second  order  for  the  expectation 
value  of  Q  = 
writing 


Q^(r^)  can  illustrated  from  Equation  (7 6)  by 
i 


<Q> 


approximate 


G>*  Q$  dT 


(77) 


=  <Q>eXact  ~  f  r*(^dT  +  f 

+  £2[<Q>exact  +  /  X^XdT) 


which  gives 


<Q>  .  ^  =  <Q>  „  +  €2[<Q>  +  /  X*QXdr]  (78) 

approximate  exact  exact  J 

because  the  integrals  <X |Q |¥>  reduce  to  a  sum  of  integrals 

®*(r....r  )Q(r.)  )dT  ...dT  (79) 

v  1  n/xv  1  v  1’  n  1  n 

each  of  which  then  vanishes  by  application  of  Brillouin's  theorem. 

It  is  generally  accepted  that  Hartree-Fock  wavefunctions 
define  a  well -averaged  environment  for  the  individual  electrons,  making 
them  suitable  for  the  determination  of  one-electron  properties.  How¬ 
ever,  it  does  not  necessarily  follow  that  all  one-electron  operators 
will  yield  expectation  values  with  the  same  degree  of  accuracy. 


<xiqii>  =  y  f 
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Although  the  one-electron  operators  only  depend  upon  an  average  dis¬ 
tribution,  a  particular  one-electron  operator  may  act  with  a  different 
weight  on  different  space  regions  over  which  the  wavefunction  is  defined. 
If  the  function  is  poorly  defined  over  the  principal  regions  on  which 
the  operator  acts,  then  the  expectation  values  will  be  poor.  The  domain 
of  defintion  of  the  operator,  \  exp  (iK.s  •  r . )  ,  giving  the  expectation 


j 

values  for  the  coherent  scattering  factors  and  intensities  is  seen  to 
be  the  same  as  that  for  the  total  wavefunctions  from  which  they  were 
computed,  namely  0  <  cp  <  a  ,  0  <  @  <  2jt,  0  <  r  <  oo  .  It  was  previously 

shown  that  the  coherent  scattering  was  directly  related  to  the  radial 

■& 

distribution  function,  a  function  wich  is  well  represented  in  the 
Hartree-Fock  approximation. 


In  view  of  Brillouin's  theorem  and  the  discussion  regarding  the 
accuracy  of  one-electron  properties  calculated  from  approximate  wave- 
functions,  it  is  evident  that  the  accuracy  of  the  coherent  scattering 
factors  listed  in  the  appendices  is  dependent  upon  the  extent  to  which 
the  analytical  functions  represent  the  true  state  functions  ¥  ,  A 
convenient  measure  of  the  deviation  from  exact  values  of  one-electron 
properties  calculated  from  the  analytical  functions  is  provided  from 
the  magnitude  of  e  .  As  seen  from  Equation  (76),  an  accurate 

determination  of  e  presupposes  a  knowledge  of  the  exact  state  function 

** 

¥  .  Estimates  of  e  should  therefore  be  made  before  passing  judgment 


For  a  descussion  of  radial  probability  plots,  reference  cajq.  be  made 
to  D.  R.  Hartree,  The  Calculation  of  Atomic  Structures  [19]. 


From  the  coefficients  of  configuration  functions,  Goodisman  and 
Klemperer  [18]  estimate  at  .045  for  He  . 
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on  the  ability  of  these  analytical  functions  to  yield  accurate  one- 
electron  expectation  values.  In  addition,  a  comparison  of  several  one- 
electron  properties  computed  from  these  functions  should  be  made  with 
good  experimental  values. 

Until  accurate  experimental  values  for  the  coherent  scattering 
factors  become  available  or  until  exact  wavefunctions  for  all  the  atomic 
systems  considered  in  the  present  computations  are  known,  no  general 
attempt  can  be  made  to  assess  the  goodness  of  the  present  results.  How¬ 
ever,  the  present  values  for  the  coherent  quantities  are  as  good  as 

■%r 

existing  values  computed  from  numerical  Hartree-Fock  functions .  The 
results  for  the  coherent  quantities  tabulated  in  the  appendices 
represent  Hartree-Fock  expectation  values  accurate  to  at  least  four 
significant  figures  (with  the  fifth  speculative). 

The  total  scattered  intensity  was  seen  to  depend  on  the 
distribution  function  of  one  electron  with  respect  to  another,  that  is, 
on  electron  correlation.  As  the  analytical  Hartree-Fock  functions  do 
not  account  for  electron  correlation,  any  two-electron  properties 
computed  from  them  would  be  expected  to  be  in  error.  The  extent  of 
the  error  associated  with  the  incoherent  scattering  factor  will  depend 
on  the  relative  weight  of  the  contributions  of  the  two-electron  and 
one-electron  operators  to  the  expectation  value.  This  is  readily 


Rather  extensive  tabulations  are  required  to  demonstrate  this  point. 
Comparison  of  present  results  can  in  some  cases  be  made  with  those 
values  computed  from  numerical  Hartree-Fock  functions.  In  particular, 
the  papers  of  A.  J.  Freeman  [15,16]  provide  results  for  some  of  the 
metals  and  metallic  ions  of  the  first  row  transition  elements. 
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illustrated  by  considering  the  equations  for  the  total  scattered 
intensity. 


i/ici  =  An)  +  s(n) 


sin  dt  -i-  N  }  , 

lit  J 


(80) 


and  the  equation  for  the  coherent  intensity, 


sin^ 

|ir 


(81) 


Solving  for  S(|i)  gives 

pOO  -00 

S(n)  =  Ii/Ici  =  N  +  /  p(t)  sin  dt  '  [  /  D(r)  sin  fjjj  dr]2*  (82 

Equation  (82)  states  that  the  expectation  value  for  S(|i)  is  given  as 
the  difference  of  the  total  intensity,  a  two-electron  property,  and  the 
square  of  a  one-electron  property,  the  coherent  scattering  factor,,  How¬ 
ever,  as  the  number  of  electrons  increases,  the  coherent  intensity  rapidly 
approaches  the  total  intensity,  thereby  implying  that  the  dominant 
contributions  to  the  incoherent  scattering  factor  (intensity)  are  due 
to  the  one-electron  operator.  The  dominant  contribution  to  the  incoherent 
intensity  comes  from  the  two-electron  operator  only  if  the  scattering 

system  contains  a  relatively  small  number  of  electrons.  The  following 

2  2 

table  lists  the  ratio,  [F  (p.)  -{-  S(|a)]/F  (p)  of  the  total  to  the 
coherently  scattered  x-ray  intensities  for  the  neutral  atoms  He  -  Kr  , 
calculated  at  a  value  of  (sin  a) =  0.300  „ 


nun  aOt  t«  mf 

aub  3-5  sXa\  .:;e-jD  r,  30  ool  sftistrj;  oe  tJnais  Joonl  or, 3  oJ  afloiir  ,i*-  f0o 
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Table  2.  Ratio  of  Total  to  Coherent  Scattered  X-Ray  Intensities 


for 

the  Atoms 

He  - 

Kr 

Atom 

Ratio 

At  om 

Ratio 

Atom 

Ratio 

He 

2.2793 

Si 

I.O856 

Fe 

1 .0298 

Li 

1.7992 

P 

1.0915 

Co 

1.0268 

Be 

1.8712 

S 

1.0828 

M 

1.0240 

B 

1.8223 

Cl 

1.0777 

Cu 

1.0217 

C 

1.60032 

Ar 

1.0709 

Zn 

I.OI96 

N 

1.3870 

K 

1.0637 

Ga 

I.OI87 

0 

I.2532 

Ca 

1.0584 

Ge 

1.0182 

F 

I.I67I 

Sc 

I.O527 

As 

I.OI79 

Be 

1.1156 

Ti 

1.0472 

Se 

I.OI77 

Na 

1.0934 

V 

1.0420 

Br 

1.0174 

Mg 

1.0861 

Cr 

1.0374 

Kr 

1.0171 

A1 

1.0847 

Mn 

1.0332 

At  least  for  the  heavier  elements,  Table  2  clearly  indicates 
that  the  major  portion  of  the  total  scattered  intensity  corresponds  to 
radiation  which  has  been  elastically  scattered,.  For  the  case  of  oxygen, 
the  ratio  of  the  cement  to  the  total  intensity  is  already  80f0  ,  for 
the  present  value  of  the  scattering  variable.  At  a  smaller  scattering 
angle,  the  contribution  of  the  coherent  intensity  to  the  total 
intensity  would  be  even  greater^  and  in  the  limiting  case,  (sin  a)/A  =  0  , 
the  total  intensity  would  be  completely  represented  by  the  coherent 


contribution . 


A 

-  J,  ■ 
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It  is  noticed  that  in  going  from  Li  to  Be  ,  the  ratio 
increases,  as  it  does  for  the  sequence  A1  -  Si  -  P  .  The  general 
trend  is  to  a  decreasing  ratio  with  increasing  number  of  electrons, 
implying  that  coherent  contributions  to  the  total  intensity  are  indeed 
dominant.  These  apparent  anomalies  must  be  due  to  an  increase  in  the 
contribution  of  the  incoherent  scattering  to  the  total,  but  in  view 
of  the  fact  that  electrons  are  being  added  to  partially  filled  shells, 
the  behaviour  may  possibly  be  interpreted  as  a  manifestation  of  electron 
correlation*  Similarily,  the  change  in  the  ratio  on  going  from  Ne  to 
Na  is  *0222  ,  but  is  only  *0073  for  Na  to  Mg  *  A  lesser  effect 
is  observed  for  the  sequence  Ar  -  K  ,  K  -  Ca  ,  where  the  ratio  changes 
are  *0072  and  .0053  >  respectively.  Similar  trends  are  noticed  in 
the  changes  Ca  -  Sc  ,  Sc  -  Ti  ,  and  Zn  -  Ga  ,  Ga  -  Ge  „  In  all 
cases  the  change  in  the  ratio  is  more  pronounced  when  a  shell  begins 
to  be  occupied  with  the  change  then  being  slower  as  the  occupation  of 
the  shell  increases.  The  same  effect  is  noted  in  the  3d  shell  of 
the  transition  elements  where  the  change  in  ratio  monotonically  decreases 
until  the  shell  is  half-filled  at  Mn  ,  increases  from  Mn  to  Fe  , 
and  then  monotonically  decreases  until  the  shell  is  filled  at  Zn  . 

In  the  absence  of  the  electron  correlation,  the  ratio 
differences  should  form  a  uniformly  decreasing  sequence  from  He  to 
Kr  .  That  fluctuations  within  a  given  shell  are  not  as  marked  as  those 
encountered  on  going  from  a  filled  shell  to  an  open  shell  with  one 
electron  is  taken  as  evidence  for  the  effects  of  electron  correlation 


on  the  total  scattering# 
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The  innermost  electrons,  especially  in  the  heavier  elements, 
contribute  most  heavily  to  the  coherent  and  hence  the  total  scattered 
intensity.  On  the  other  hand,  the  increases  in  the  ratios  for  Li  -  Be 
and  A1  -  Sc  -  P  serve  to  illustrate  that  contributions  from  outermost 
shells  are  important  to  the  incoherent  intensity. 

The  ratios  listed  in  Table  2  indicate  that  incoherently 
scattered  radiation  corresponds  to  a  significant  portion  of  the  total 
scattered  radiation,  at  least  for  the  lighter  atoms.  Although  the  ratios 
were  calculated  with  values  determined  from  approximate  wave func t ions , 
Table  2  gives  an  indication  of  the  magnitude  of  the  necessary  corrections 
to  be  applied  to  coherent  scattering  factors  determined  from  experimental 
measurements  of  total  scattered  intensities,,  In  order  to  assess  the 
goodness  of  the  present  results  for  the  incoherent  scattering  factors, 
comparison  will  be  made  with  the  expectation  values  calculated  by 
Bartell  and  Gavin  [1,  2]  from  a  correlated  function  which  accounts  for 
92.1$>  of  the  correlation  energy.  The  comparison  will  be  made  for 
helium  only,  since  the  above  table  indicates  that  for  this  atom  the 
incoherently  scattered  radiation  contributes  the  major  portion  to  the 


total. 


- 

' 


■ 
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Table  5-  Incoherent  Scattering  factors,  S(|a)  ,  for  Helium 


0  _  1 

1  a)/A,  A 

a 

Hartree-Fock 

tP(r1)'P(r2) 

Roothaan-Weiss^ 

cp(ri)cp(r2)x(r12) 

Relative  Err 

0.025 

O.O21579 

0.02052 

5.15 

0.050 

0.084754 

0.08032 

5.39 

0.075 

0.18421 

0.1744 

5.61 

0.100 

0.31224 

0.2954 

5.70 

0.125 

0.45981 

0.4348 

5-75 

0.150 

O.6I780 

0 . 5840 

5.80 

0.175 

O.77797 

0.7355 

5.74 

0.200 

0.93363 

0 . 8830 

5.72 

0.250 

1.2136 

1.1485 

5.68 

0.300 

1.4380 

1 . 3669 

5.20 

0.400 

1 .7274 

1.6621 

3.92 

0.500 

1.8705 

1.8219 

2.67 

0.600 

1.9377 

I.9052 

1.70 

0.700 

1.9692 

1.9485 

1.06 

0.800 

1 . 9843 

1.9712 

0.665 

Present  calculation,  using  the  analytical  Hartree-Fock  function 
(0$  of  correlation  energy)  of  Clementi  [8], 

92.1$>  of  correlation  energy,  reference  [42], 


The  results  of  the  present  calculation  are  in  relatively  good 


agreement  with  those  from  the  correlated  function,  showing  a  maximum 


• 

, 

■ 
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relative  error  of  5.80 $  at  (sin  a)/A  =  0.175* 

For  helium,  it  is  observed  that  the  values  of  S(|_i)  computed 
from  the  analytical  function  are  larger  than  those  evaluated  from  the 
correlated  function.  However,  any  statement  regarding  the  accuracy  of 
the  values  of  S(p.)  for  the  remaining  atoms  Li  -  Kr  must  again  be 
based  on  a  comparison  with  values  computed  from  accurate  functions, 
or  on  a  comparision  with  experimental  results.  As  neither  experimental 
values  nor  exact  wavefunctions  are  available  for  all  the  systems  considered 
in  these  calculations,  a  general  assessment  of  the  goodness  of  the  values 
of  S(n)  as  calculated  from  dementi’s  functions  is  not  possible. 

The  discussion  of  the  results  for  the  incoherent  terms  involved 
values  computed  from  Equation  (25),  which  is  in  effect  a  zero-order 
equation  with  respect  to  contributions  from  excited  states.  Reference 
can  be  made  to  the  work  of  Bonham  [5]  for  corrections  to  third  order. 

As  pointed  out  by  Bonham,  these  corrections  are  important  in  the  case 
of  low-energy  incident  radiation  or  heavy  atoms,  where  the  first 
ionization  potentials  are  low.  The  effects  of  chemical  bonding 
on  the  scattering  factors  has  been  discussed  by  McWeeny  [37]° 

Polarization  and  relativistic  effects  have  not  as  yet  been  completely 
assessed . 

Even  though  Equations  (24)  and  (25)  of  the  theory  section 
were  derived  for  gaseous  atoms,  they  can  be  applied  to  crystals  and 
molecules  if  the  appropriate  many-electron  state  functions  for  these 
systems  are  known.  Until  these  functions  become  available,  the  free 


J  :  , 
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atom  scattering  factors  will  have  to  be  used  when  dealing  with  molecules 
and  crystals.  For  crystals  with  large  numbers  of  well-ordered  atoms, 
the  coherent  scattering  factors  given  here  should  be  of  particular  value 
in  the  calculation  of  total  intensities.  Correlation  effects  in  crystals 
would  not  be  expected  to  extend  beyond  nearest  neighbors,  and  for  well- 
ordered  crystals  their  effects  should  be  insignificant. 
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Explanation  of  Appendices 

Table  1  of  Appendix  1  collects  the  values  of  the  coherent 
electron  scattering  factor,  f^^(O),  for  the  positive  ions  neutral 
atoms  and  negative  ions  He  to  Kr  .  Table  2  of  Appendix  1  lists  the 

values  for  the  coherent  x-ray  scattering  factor,  F  ,  coherent 

2 

x-ray  scattered  intensity,  F  ,  incoherent  electron  and  x-ray  scatter¬ 
ing  factor  S,  and  the  square  of  the  electron  scattering  factor, 

(N  -  F),  for  all  the  neutral  atoms  He  to  Kr  .  The  column  at  the 

left  hand  of  each  page  lists  the  range  of  values  of  the  scattering 

° 

variable,  (sin  <x)/A  ,  in  A  )  for  which  the  evaluations  were 
performed . 


The  remaining  tables  of  Appendices  2  and  3  collect  values 
for  the  coherent  and  incoherent  scattered  x-ray  intensities  [a],  The 
column  headings  are  as  follows: 

2 

Coherent  scattered  x-ray  intensity'  F  =  (SF(l,l))  **  2 
Incoherent  scattered  x-ray  intensity,  S  =  Z  -  F  [b] 

Sum  of  the  squares  of  the  diagnol  elements®  S(F(l,l))  **  2  [c] 

Sum  of  the  squares  of  diagnol  and  non-diagnol  elements;  2  =  S(F(l,j))  **  2 
Ratio  of  the  incoherent  to  coherent  x-ray  scattered  intensity;  INC/COH 

The  tables  of  Appendix  2  list  the  above  values  for  positive 
ions,  whereas  the  tables  of  Appendix  3  present  values  for  negative  ions. 

The  range  of  the  scattering  variable,  (sin  a)/A  ,  is  the  same  for  all 
tables  of  Appendices  2  and  3  as  that  for  Table  2  of  Appendix  1, 
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[a]  The  corresponding  values  for  the  coherent  x-ray  scattering  factors 

and  the  coherent  electron  scattering  factors  and  intensities  can  be  found 

2 

from  the  values  of  F 

[b]  This  F  term  corresponds  to  the  complete  summation  term,  Z  , 
of  Equation  (25) . 

2 

[c]  This  is  the  summation  term,  Z^  |fjj  |  ,  of  Equation  (25). 
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Appendix  1 


Table  1. 


Atomic  factor  '(0)  (in  A  )  electron  scattering  for 

a) 

positive  ions,  neutral  atoms,  and  negative  ions. 


Postive 

Ion 

Neutral 

Atom 

Negative  Ion 

stem 

State 

f^L)(0) 

State 

State 

f7}(°) 

He 

xs 

0.42 

Li 

xs 

0,16 

2s 

3.29 

xs 

15.51 

Be 

2s 

1.16 

1s 

3.05 

B 

ls 

1,40 

2P 

2.79 

3p 

7.02 

C 

2P 

1.42 

5p 

2.43 

ks 

4.64 

N 

5P 

1.40 

2.13 

3P 

3.93 

0 

2s 

1.30 

5P 

1.97 

2'P 

3.30 

F 

3P 

1.24 

2P 

1.81 

xs 

2.81 

Ne 

2 

P 

1.20 

xs 

1.65 

Na 

xs 

1.13 

2s 

4.79 

h 

17.42 

Mg 

2s 

2.50 

5.22 

A1 

ls 

3.07 

2P 

5.90 

3p 

13.08 

Si 

2P 

3.60 

3.p 

5.69 

4S 

9.60 

P 

5P 

3.70 

5.34 

5P 

8.58 

S 

4s 

3.70 

5P 

5.14 

2P 

7.58 

Cl 

3p 

3.64 

2P 

4.87 

ls 

6.75 

A 

2 

P 

3.60 

*8 

4.59 

K 

xs 

3.44 

2 

S 

9.03 

h 

28.14 

• 

Table  1  (continued) 


Positive  Ion  Neutral  Atom  Negative  Ion 


System 

State 

f(1)(0) 
e  x 

State 

f^(0) 

State 

f^CO) 

Ca 

2s 

5.51 

Ls 

9.98 

Sc 

5F 

4.04 

2d 

9.38 

^F 

19.63 

Ti 

V 

3.92 

5p 

8.85 

4f 

19.19 

V 

3.81 

8.39 

5d 

18.31 

Cr 

6S 

3.70 

7S 

6.70 

6s 

17.50 

Mn 

5D 

3.63 

6s 

7.60 

5d 

16.28 

Fe 

s? 

3.54 

5d 

7.26 

^F 

15.56 

Co 

5F 

3.44 

kF 

6.94 

5f 

14.63 

Ni 

2d 

3.40 

3f 

6.67 

2s 

8.57 

Cu 

Ls 

3.30 

2s 

5.90 

7s 

13.38 

Zn 

2s 

4.00 

xs 

6.20 

Ga 

Ls 

4.50 

2P 

7.21 

3P 

14.78 

Ge 

2P 

5.10 

3P 

7.33 

11.42 

As 

3P 

5.49 

7.24 

3p 

10.80 

Se 

4s 

5.50 

3P 

7.30 

2P 

10.09 

Br 

5P 

5.62 

2P 

7.14 

9.39 

Kr 

2P 

5.64 

Ls 

7.00 

Values  calculated  from  the  analytical  Hartree-Fock  functions  of 
Clementi  [8]  and  Malli  [  36]  .  Strand  and  Bonham  [46]  have  also 
calculated  these  values  for  neutral  atoms. 


Appendix  1  Table  2 


HE  IS  (2)  IS 


F 

S 

F**2 

(N-F  )** 

2 

0.000 

0.2C000E 

01 

-0.93937E- 

04 

0.40002E  01 

0. 22060E 

-08 

0.  025 

0. 19892E 

0  1 

0. 21 579E- 

01 

0.39568E  01 

0.  1  1705E 

-03 

0.050 

0. 19572E 

01 

C. 84754E- 

01 

0.3830 5 E  01 

0.  18349E 

-02 

0.075 

0. 19057E 

01 

0.  184 2  IE 

CO 

0.36316E  01 

0. 88982E 

-02 

0.  100 

0,  1  8 3  73 E 

01 

0 . 3 l  224  E 

CO 

0.33755E  01 

0.264 84 E 

-01 

0.  125 

0. 17551E 

0  1 

0. 4  598  1  E 

CO 

0.30804E  01 

0. 59975E 

-01 

0.  150 

0. 16626E 

01 

0. 61 78CE 

CO 

0.276 4 4E  01 

0. 11381E 

00 

0.  175 

0. 15633E 

01 

C. 77797E 

CO 

0.24441E  01 

0. 19067E 

00 

0.200 

0. 14 60  4 E 

01 

0 . 9  3  3  6  3  E 

CO 

0.21327E  01 

0.291  1 8  E 

0  0 

0.  225 

0. 13566E 

0  1 

0. 10799E 

Cl 

0 . 1  84 03 E  01 

0.41402E 

00 

0.250 

0. 12  54  IE 

01 

0. 12136E 

01 

0.15729E  01 

0. 55629E 

00 

0.27  5 

0.1 1  54 9 E 

01 

0  .  1  3  3  3  1 E 

01 

0.13339E  01 

0.714  14E 

00 

0.  300 

0. 1 0602  E 

01 

0. 1 4330E 

01 

0.1  1 240  E  01 

0. 88323E 

00 

0.  325 

0 .97080E 

00 

0. 1 52  88  E 

01 

0.94246E  00 

0.  1 05  92  E 

0  1 

0.350 

0.68725E 

00 

0. 16064E 

01 

0.78722E  00 

0. 12382E 

01 

0.375 

0.8C976E 

00 

0. 16721E 

01 

0. 6  55  72  E  00 

0. 14167E 

0  1 

0.4  00 

0. 7  3  834  F. 

00 

6.1  7 2  74  E 

C  1 

0.54515E  00 

0. 15918E 

01 

0.4?5 

0 . 67  28  4  E 

00 

0. 17736E 

01 

0.45271E  00 

0.  1  7 6 1 4  E 

01 

0.4  50 

0.61 298E 

00 

0. 18121E 

Cl 

0.37575E  00 

0.  192  38E 

0  l 

0.475 

0. 55H46E 

CO 

0 .  1  8  4  4  1 E 

01 

0.31188E  00 

0. 2U780E 

01 

0.500 

C.50891E 

00 

0. 187G5E 

01 

C.25899E  00 

0 . 222  33 E 

01 

0. 550 

0 . 4  2  320  E 

00 

0. 19104E 

01 

C.17910E  00 

0. 24863E 

01 

0.600 

0.35239E 

00 

0. 19377E 

01 

0.12453E  00 

0.271 30E 

01 

0.650 

0. 29528E 

00 

0. 19564E 

01 

0.871 92 E-01 

0. 29061E 

0  1 

0.  700 

0.2480 4 E 

00 

0. 19692E 

01 

0.61525E-01 

0. 30694E 

01 

0.750 

0. 20923E 

00 

0.  19  78  IE 

01 

0.43776E-01 

0. 32069E 

01 

0.800 

0. 17  724E 

00 

0. 19843E 

01 

0.3 14 15E-01 

0. 33224E 

01 

0.850 

0. 15C80E 

00 

0. 19886E 

01 

0.22740E-01 

0. 34196E 

01 

0.900 

0. 12885E 

00 

0. 1991  YE 

01 

0. 16602E-01 

0 . 3  50  1 2  E 

01 

0.950 

0. 1 1056E 

00 

0. 19939E 

01 

0. 12223E-CI 

0. 35700E 

01 

1.000 

0.95255E 

-01 

0. 19955E 

01 

0.90736E-02 

0. 36281 E 

0  L 

1.  IOC 

0.71555E 

-01 

0. 19974E 

01 

0.5 1201E-02 

0.37189E 

01 

1.200 

0.54562E 

-.01 

0. 19985E 

01 

0.29770E-02 

0. 37847E 

01 

1.300 

0.42  136E 

-01 

0. 1 999 1 E 

01 

0.1 7797E-02 

0.38330E 

01 

1.400 

0.33040E 

-01 

0. 19995E 

01 

0. 10916E-02 

0.38689E 

01 

1.500 

0.26186E 

-01 

0. 19997E 

01 

0.68569E-03 

0. 38959E 

01 

1.600 

0.20982E 

-01 

0. 19998E 

01 

0.44025E-03 

0 . 39 1 65  E 

01 

1.700 

0 . 1 6  9  8  4  E 

-01 

0. 19999E 

01 

0.28845E-03 

0. 39324E 

0  l 

1.8C0 

0.13877E 

-01 

0. 19999E 

01 

0. 19256E-03 

0. 39447E 

01 

1.9CC 

0. 1 1  437 E 

-01 

0. 19999E 

01 

0.1 3080E-03 

0.39544E 

01 

2.000 

0.95020E 

-02 

0.20 COO E 

01 

0. 90288 E -04 

0. 39621E 

01 

2.  ICO 

0 .79536E 

-02 

0. 2CC00E 

01 

0.63260E-04 

0.39682E 

01 

2.200 

0.6703 9 E 

-02 

0. 2CCC0E 

01 

0.44942E-04 

0 . 397  32  E 

0  1 

2.  300 

0.56872E 

-02 

C.20C00E 

01 

0. 3234 4 E -04 

0.397  73 E 

0  1 

2.400 

0.48538E 

-02 

0.20G00E 

01 

0.23560E-04 

0. 39806E 

0  1 

2.50C 

0.4  1660 E 

-02 

0. 2C0CCE 

01 

0. 1 7356E-04 

0.  398  34  E 

0  l 

2.700 

0.31  1 69  E 

-02 

0 . 2  0  0  0  0  E 

01 

0.971 52E-05 

0.39875E 

01 

2.900 

0. 23758E 

-02 

0.2000CE 

Cl 

0.56444E-05 

0. 39905E 

0  l 

3.  100 

0. 184  10E 

-02 

0 . 20C00E 

01 

0.33891E-C5 

0. 39926E 

01 

3.  300 

0 . 14  476E 

-02 

0. 20000E 

01 

0 . 20956E -05 

0.399426 

0  1 

3.500 

0.1  1  534E 

-02 

0. 20CC0E 

01 

0. 13302E-05 

0. 39954E 

01 

■ 

* 

. 
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LI  IS  (2 ) 2S  (1)  2S 


F 

s 

F**2 

( N-F ) ** 

2 

0.000 

0. 30000E 

01 

-0.41574b- 

04 

0.90001b 

01 

0. 43T49b 

-09 

0.025 

0. 29177E 

01 

0.  15851b 

00 

0.85129b 

01 

0.67763b 

-02  N 

0.050 

0.27075b 

01 

0. 50856b 

00 

0. 73305b 

01 

0.85562b 

-01 

0.075 

0 . 2  4  50  l  L 

01 

0. 83532E 

00 

0. 60028E 

01 

0.30244b 

00 

0.100 

0.221 50E 

01 

0.  104  76b 

01 

0.49060b 

01 

0.61629b 

00  ' 

0. 125 

0. 2033  lb 

01 

0.11737b 

01 

0.41336b 

01 

0.93433b 

00 

0.150 

0.  19036b 

01 

0. 12651b 

01 

0. 36237b 

01 

0.  1  2  02  1  L 

01 

0.175 

0.18116b 

01 

0.13512b 

01 

0. 32819E 

01 

0.14123b 

01 

0.200 

0.  1 7  4 1 5  E 

01 

0. 144  12E 

01 

0.30330b 

01 

0. 15837b 

01 

0.225 

0.  16820b 

01 

0.15358b 

01 

0.28290b 

01 

0.17372b 

0  1 

0.2  50 

C. 16260E 

01 

0. 16331b 

01 

0.26437b 

01 

0. 188805 

01 

0.2  75 

0.  15700b 

01 

0. 1731 3E 

01 

0. 24648E- 

01 

0.20450b 

01 

0.300 

0. 15126b 

01 

0. 18286E 

01 

0.22879b 

01 

0.22124b 

0  l 

0.325 

0.14535b 

01 

0.19237b 

01 

0.21127b 

01 

0.23917c 

0  1 

0.350 

0. 13931b 

01 

0. 201 36b 

01 

0.19406b 

01 

0. 25823c 

01 

0.37  5 

0.  13318E 

01 

0.210  .b 

01 

0. 17736E 

01 

0.27830b 

01 

0. 400 

0.  12703E 

01 

0.218 

01 

0.  16136b 

01 

0.29919b 

01 

0 .  A 25 

0.12091b 

01 

0.226 

01 

0. 14620E 

01 

0.32072b 

01 

0.4  50 

0.  11488b 

01 

0.233V  lb 

01 

0. 13199E 

01 

0. 34268E 

01 

0.475 

0.  10898E 

01 

0.24042b 

01 

0. 11877E 

01 

0.36487b 

0  1 

0.500 

0.  10324E 

01 

0. 24638b 

01 

0. 10659E 

01 

0.38713b 

0  l 

0.550 

0.92345b 

00 

0.25731b 

01 

0. 85277b 

00 

0.43120b 

01 

0.600 

0. 82319b 

00 

0. 26610b 

01 

0. 67764b 

00 

0. 4 /385b 

0  l 

0.650 

0. 73215E 

00 

0.27319E 

01 

0. 53604 E 

00 

0.51432b 

0  1 

0.700 

0. 65028b 

00 

0.27885b 

01 

0.42286b 

00 

0.55212c 

01 

0.750 

0.57719b 

00 

0.28334b 

01 

0.33315b 

00 

0.58700b 

01 

0.800 

0. 51227E 

00 

0. 28688E 

01 

0.26242E 

00 

0. 61880b 

0  l 

0.850 

0.45483E 

00 

0.28966b 

01 

0. 20687b 

00 

0. 64779E 

0  1 

0.900 

0.40413E 

00 

0.29183b 

01 

0. 16332b 

0  0 

0.67385b 

0  1 

0.950 

0.359 4 5E 

00 

0.29354b 

01 

0. 12921 t 

00 

0.69725b 

01 

1.000 

0.3201  IE 

00 

0.29488b 

01 

0. 10247b 

00 

0. 71B18E 

01 

1.  100 

0. 25499E 

00 

0.29675b 

01 

0.65020E 

-01 

0.75351b 

01 

1.200 

0.204 4 4E 

00 

0.29791b 

01 

0.41797b 

-01 

0.78151b 

01 

1.300 

0.16506b 

00 

0.29864b 

01 

0. 27245b 

-01 

0.80369b 

01 

1.400 

0.  1 342  1 E 

00 

0.29910E 

01 

0. 18013b 

-01 

0. 82127b 

01 

1.500 

0. 1G990E 

00 

0.29940b 

01 

0. 12079E 

-01 

0.83527b 

0  l 

1.600 

0.90621E 

-01 

0. 29959E 

01 

0.82122b 

-02 

0.84645b 

01 

1.700 

0.75225b 

-01 

0.29972b 

01 

0. 56588b 

-02 

0.85543b 

01 

1.800 

0. 62848E 

-01 

0. 29980E 

01 

0.39499b 

-02 

0.862695 

01 

1.900 

0. 52G32E 

-01 

0.29986E 

01 

0.27912E 

-02 

0. 86658E 

01  ; 

2.000 

0 . 4  4  6  7  5  E 

-01 

0.29990b 

01 

0. 19958E 

-02 

0. 87339E 

01 

2.100 

0.37909b 

-01 

0. 29993b 

01 

0.14431b 

-02 

0. 87735b 

01 

2. 200 

0.32475E 

-01 

0.299 9 5E 

01 

0.10546b 

-02 

0. 880d2  b 

01 

2.300 

0. 27903E 

-01 

0. 29996b 

01 

0.77857b 

-03 

0.88334b 

0  1 

2.400 

0.24089E 

-01 

0.29997b 

01 

0.58029b 

-03 

0. 88560E 

01 

2.500 

0  •  20c-92b 

-01 

0. 29998b 

01 

0.43646b 

-0  3 

0.8875  lb 

01 

"2.700 

0.  15915b 

-01 

0. 29999b 

01 

0. 25328E 

-03 

0. 8 904 8 E 

01 

2.900 

0. 12314E 

-0  1 

0.29999b 

01 

0. 15162E 

-03 

0.69263b 

0  l 

3.100 

0. 96624 E 

-02 

0. 30000E 

01 

0.93362b 

-04 

0.89421b 

0  l 

3.300 

0.76792b 

-02 

0.30000b 

01 

0.58970b 

-04 

0.89540b 

01 

3.500 

0.61  7 42E 

-02 

0. 30000b 

01 

0. 38 120b 

-04 

0. 89630c 

01 

_  J 


I 
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BE  IS  (2 ) 2S  (2 ) 


IS 


F 

s 

F**2 

(  M  —  F  )  £  ■>' 

2 

•  0.000 

0.4000  1  E 

01 

-0. 13191E- 

0  3 

0. 1600  IE 

02 

0.43450b 

-08 

.  0.025 

0.39220E 

01 

0. 15316E 

CO 

o: l  5  3  3  2  E 

02 

0. 60073E 

-02 

0.050 

0 . 3  7  067E 

01 

0. 54602E 

CO 

0. 1  3  7  39  E 

02 

0. 86038E 

-01 

0.075 

0. 34021E 

01 

0.  10279E 

01 

0.1 1575E 

02 

0. 35744E 

00 

0.  100 

0.30 65  IE 

01 

0. 1462SE 

01 

0.93948E 

01 

0. 87405E 

GO 

0.  125 

0 . 2  7  4  2  4  E 

0  1 

0.  1 7  865  E 

01 

0. 75209E 

01 

0.  15815b 

0  l 

0.150 

0.24624E 

01 

C. 20C04E 

01 

0. 60635E 

0  1 

0. 23642E 

0  l 

0.175 

0.2235 4 E 

01 

0. 21 364 E 

01 

0.49971 E 

01 

0.311 37b 

01 

0.200 

0. 20596E 

01 

C. 22281E 

01 

0 . 4242 1 E 

01 

0. 37651 E 

01 

0.225 

0. 19269E 

01 

0.22995E 

01 

0 . 3  7  1  3  1  E 

01 

0.42976E 

01 

0.250 

0. 18274E 

01 

0. 2 3643 E 

01 

0. 33395E 

01 

0.47201b 

01 

0.275 

0. 17517E 

01 

0.24284E 

01 

0. 30685E 

01 

0.  5 0 5 4 8 E 

01 

0.300 

0. 1692  IE 

01 

0. 24942E 

Cl 

0. 2  8632  E 

01 

0. 53264E 

01 

0.  325 

0. 1642  7E 

01 

0 . 2  5  6  1  6  E 

Cl 

0.26986E 

01 

0.  5 5 5 6 7 E 

0  1 

0.350 

0 . 15994E 

01 

0. 26304E 

01 

0.25581E 

01 

0. 57629E 

01 

0.37  5 

0. 15592E 

01 

0. 26998E 

01 

0.24312b 

01 

0.59574b 

01 

0.400 

0. 15204E 

01 

0.27692E 

01 

U.231  15E 

01 

0.61486b 

01 

0  .  4  2  5 

0 . 1 4  8  1  7  E 

01 

0.28380E 

01 

0.2  1954 E 

01 

0.63419E 

01 

0.450 

0. 14426E 

01 

0. 29058E 

01 

0 . 2  0  8  1  0  E 

01 

0 «  6  5  4  0  4  E 

01 

0.475 

0.  l*i023E 

01 

0. 29722E 

0  1 
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